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C/!! \ Abstract 

Theory of diffusion-mediated reactions is aheady estabhshed for the target problem in the dilute 
limit, where the immobile target is surrounded by many quenchers. For lattice random walks in 
^ • the crowded situation, each quencher is surrounded by other quenchers differently. As a result, each 



I 

o 



quencher migrates differently in the presence of site blocking effects. However, in the conventional 
theory, such difference is ignored and quenchers are assumed to move independently of each other. 
In this paper, theory of diffusion-mediated reactions of target problem is developed by taking into 
account the site blocking effects for quencher migration and the difference in the configuration of 
quenchers around each quencher. Our result interpolates between those in high and low limits of 
quencher concentrations and is a lower bound of the survival probability. In the static limit, the exact 
result is reproduced for a localized sink. In the presence of diffusion, the approximation is better when 



(N 

, intrinsic reaction rates are low. 

!>■ 
00 
(N 
(N 

00 

O 

0\ 

o 



X 



* Electronic address: |k-seki@aist.go.jp 



Electronic address: . t achiy a@ aist . go . j p| 



1 



I. INTRODUCTION 



Theory of diffusion-mediated reactions is already established for the target problem in the 
dilute limit. Consider reaction between excited probe A* and quencher B which 

deactivates excited probe. In ordinary experimental conditions, A* is minority species and B 
is majority species. The case in which minority species are immobile and majority species are 
mobile is in general called the target problem. The opposite case is called the trapping problem. 
In this paper, we consider the target problem. We are interested in the decay of fraction of 
A* after pulsed excitation or the fraction of A* when A is continuously excited. In the dilute 
limit, movement of each quencher can be regarded as an independent event. When the target 
is surrounded by quenchers on the lattice with M sites the survival probability of the target 
at time t denoted by PN{t) is expressed in terms of the pair survival probability at time t of a 



quencher starting from r^, f{ff,t), by ^\ 

( 



1 V 

^ exp|^-c^[l-/(f,,t)]j, (1) 

where the infinite limit of all lattice points, M, is taken and the concentration is given by, 
c = \imM~*oo N/M in the thermodynamic limit. The theory is applicable even under the long- 
range reactions and the presence of electrostatic potential among reactants, which can be taken 
into account in f{f£,t) as long as quencher concentration is dilute. 

In the conventional theory, the decay of the survival probability by bulk reactions has been 
formulated in terms of the pair probability of the target and a quencher by ignoring the ex- 
cluded volume interactions among quenchers. However, in the crowded situation, each quencher 
is surrounded by other quenchers differently. As a result, each quencher migrates differently 
in the presence of site blocking effects. In the conventional theory, such difference is ignored 
and quenchers are assumed to move independently of each other. In this paper, we take into 



account the site blocking effects for quencher migration and the difference in the configuration 
of quenchers around each quencher by applying Nakazato-Kitahara's theory of tracer diffusion. 
6| Excluded volume interactions are taken into account by prohibiting double occupancy of 
quenchers in site blocking effects. Quenchers can jump only to the empty neighbor sites. By 
noticing the success of Nakazato-Kitahara's interpolating formula on the tracer diffusion co- 
efficient between low and high concentrations of diffusing particles, we apply it to the target 
reaction on a lattice. Though site blocking is an aspect of many body interactions, rigorous 
results can be obtained by this method. In the continuous space, excluded volume interac- 
tions can be taken into account by introducing short range repulsive potentials and applying 
decoupling approximation of density correlations, as shown by Kuzovkov et al. 

II. STATIC QUENCHING WITHOUT DOUBLE OCCUPANCY 



The simplest results which prohibit double occupancy of the same site are obtained in the 
absence of diffusion. Although the results are known, we rederive them to illustrate our method. 

We consider the lattice with M sites. There are quenchers on the lattice. We indicate the 
configuration of quenchers by the set of vectors denoting the lattice sites occupied by quenchers, 
{'^i,r2, ■ ■ ■ ,r/v)- If quenchers are initially randomly distributed, the probability of taking an 
initial configuration {ri,r2, ■ ■ ■ , "Ov) is given by. 



P(ri,f2,--- ,r^;0) = 1/ (mCn) ■ 

The survival probability is obtained by applying the Cauchy's integral theorem, 8| 

If 11^ 
^^W = Tr- dx^^-—l[{l + xexp{-k{r,)t)), 

ZTTl J X M'^N „ , 



(2) 



(3) 



where the path of integration encircles the origin on the complex plane and the right hand side 
of Eq. ([3]) represents the joint probability of independent quenching events at time t from all 
possible quencher configurations. Eq. ([3]) can be rewritten as, 

M 



dx- 



X 



mCn 



exp 



In (1 + xexp(— A; (ri) t)) 



i=i 



(4) 
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In the thermodynamic hmit, Eq. (jlj) is simphfied by introducing Stirhng formula n! 
-\/2vmexp(— n)n" and applying the steepest descent method, 

M 



PN{t) 



1 

exp 



1 + x 



,M 



x^+^ mCn 



exp 



.1=1 ^ 



-k{n)t)-l) 



M 



5^1n{l + c(exp(-fc(f,)t)-l)} 



(5) 
(6) 



where c = N/M is the concentration. Eq. (jS]) is the known expression for the static quenching 
obtained by AUinger and Blumen (AB) using a different method. 9| Our method is not simple 
but shows that Eq. is correct in the thermodynamic limit for any concentration at all times 
and thus confirms the conclusion derived from the AB method, where the occupancy probability 
at each lattice site is assumed to be c instead of random occupancy of quencher sites among 
available lattice sites. In the AB method, the number of quenchers for a finite lattice with M 
sites is not necessarily equal to since the occupancy probability of each lattice site is given 
for each realization of quencher configurations. For c = 1 we find the familiar result of, 

PNit)=exp(^-Y,kin)?j, (7) 

and its Laplace transform is given by. 



Pn{s) = 1/ (^s + Y,k{re)?j 



In the opposite limit, c ~ 0, Eq. ([6]) reduces to the well-known result, 



P7v(t) ~ exp 



M 



-5^c(l-exp(-A;(f,)t)) 



(9) 



Eq. © is the static limit of Eq. (HJ. 

When reaction takes place only at a target site r^, k (r^) = A;o5r(.,rfl; the survival probability, 
Eq. ([6]), is simplified to. 



P/v(t) = 1 — c + cexp(— /cot). 



(10) 
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and the Laplace transform, P/v (s) = dtex]){—st)PN {t), is expressed as, 

p _ g + (1 - c)fco 

These trivial results will be used to check the results obtained under the presence of correlated 
diffusion. In this model, if a target site is not occupied by a quencher, the reaction never occurs 
there. However, in the presence of diffusion, even if a target site is not initially occupied by a 
quencher, a quencher may come to it by diffusion and react with it. In other words, the survival 
probability obtained under the condition of static quenching is always higher than that in the 
presence of diffusion. 

Before closing this section, we comment on the natural decay. When the natural decay of 
the target with the time constant tq is present, we multiply P/v(t) by exp(— t/ro) and s changes 
to s + I/tq in P]sf{s). The natural decay of the quencher with the time constant Tq^ is taken 
into account by replacing k (f>) by k (f>) + I/tq^. The natural decay can be taken into account 
even under the presence of correlated diffusion in the same way. 



III. QUENCHING UNDER DIFFUSION WITH EXCLUDED VOLUME INTERAC- 
TIONS 



Quenchers perform random walk on a lattice under the condition that each site cannot be 
occupied by more than one quencher at the same time. Quenchers can jump only to the empty 
neighbor sites. The movement of quencher is influenced by the position of other quenchers 
through the site blocking effects. As a result, quencher diffusion is highly correlated at high 
concentrations. In addition to correlated diffusion, reaction takes place depending on the 
distance between the quencher and the target. Since the target is immobile, the reaction rate 
depends only on the configuration of quenchers. 

As before, we consider the lattice with M quencher sites. An excited target is located at the 
origin. There are N quenchers on the lattice. If quenchers are initially randomly distributed 
over available sites, the probability of taking an initial configuration (ri, f2, • • • , tn) is given by 



5 



rN;t), 



(12) 



Eq. ([2]). The survival probability is given by, 

P^(t)=^P(fi,f2,- 

{n} 

where the summation should be taken over all possible quencher configurations. 

The self-diffusion of correlated random walk is studied by Nakazato and Kitahara in the 
ab.e„ce of .eaet.o. U S.e b.oC.g effects o. the d«us,o„ of t«ed pa.t,de ,s ealCated. 



B, 0, 



111 ] Following them, we introduce ket vectors for all accessible sites of quenchers. The 
ket vector \f, •) denotes the occupation of site r by a quencher particle, and |r, 0) represents 
that site r is an empty site. The probability of finding a configuration (f^i, ■ ■ -"nv), at time t 



averaged over all possible initial configurations of random occupation is written as 



p(fi. 



.Ht 



k=N+i / {i} \e=i / \k=N+i 



(13) 



where the sum is taken over all possible configurations of occupied sites on the M sites. H 



is given hj H = + H^^^ where describes the diffusion of quenchers. 



111 



(14) 



-f^w = ^ r/(2d) (|f„,«)(f„,0| ■ |fv„,0)(rm, •! - |f„, •)(r;,, •! ■ \rrnA){rmA\) , 

{n,rn) 

where F is the jump frequency of quencher and the sum is taken over all nearest neighbor pairs 



of accessible lattice sites by quenchers. Transition is possible from the state \ fr, 



to the 



state iTy-^ 



which indicates that the site n must be empty to accept a quencher from an 



occupied neighboring site. Similarly, if the site n is occupied by a quencher, inverse transition 



is possible for the state, 



if at least one neighboring site is vacant, if^ describes 



the reaction from an occupied site r„ with the rate k {rn), 12 



M 



(15) 



n=l 



In order to calculate the survival probability, Eq. (fT2l) . from the configuration probability, 
it is convenient to introduce the generating function, [15 1 

M 

G{x,t) = ^^xMri,r2,---,r-W;t). (16) 



The survival probability is obtained from, [15| 
where the generating function is rewritten as, 6| 

^ M M 

n ((^^' + •!) e""* n d^'^' + v^i^'^' •)) • 



(17) 



G{x,t) 

?=1 ■ fc=i 

Eq. (|T7j) with Eq. (|T8l) generalizes Eq. ([3]) by including the effect of diffusion with excluded 
volume interactions. 

In the thermodynamic limit in which M tends to infinity with the fraction of quenchers being 
fixed, c = N/M, we can apply a saddle point method to Eq. ( JT71) as we have done to obtain 
Eq. Originally, the method is introduced by Nakazato and Kitahara for the calculation 



of tracer diffusion constant of correlated random walk. 



a] The same result as theirs can be 
obtained by a different method. 16| The results of Nakazato and Kitahara is also confirmed by 
numerical simulation in 2 and 3 dimensional systems. Q, ll[ [l3] The method is based on the 
fact that the number of diffusing quenchers is conserved. In our case, the number of quenchers 
is conserved for the quencher configurations which survive reaction, and the correlated random 
walks are performed by exactly quenchers. By applying a saddle point method, Eq. ( JT71) 
becomes, 



PN{t) 



-t/ro 



M 

Hi 



M 



exp 



^0 n 



where H = exp {-6* S) H exp {6* S), S = Eii (1^^. 01 
with tan 6'* = a/c/(1 — c). H is obtained as H = Hq + Hi, where 



ri,(p)], (19) 
f>, 0)(r£, •!), and x = tan^O, 



M 



Ho = - ^k{fn) [(1 - C)|f„, •){fn, •! + C\rn, 



(20) 



n=l 



and 



M 



Hi = -^/c(r;) a/c(1 - c) i\fn, •) {fn, (p\ + |r;, 0)(f„, •!) . 



(21) 



n=l 
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By making time differentiation of Eq. f[T^ . we obtain the time evolution equation for tlie 
survival probability, 

M 



-rrPNit) = --PNit) - CKPNit) - v/c(l-c) VA: (f,) q (f,-, t) , (22) 
<it To 

where the sum of reaction rates is defined as, 

M 

K = J2k{rj)- (23) 

q {fj, t) is given by, 

q (r-, t) ^ e-*/-" (j[{re, 0| (f,-, .| j exp (i^t) |^ J] |f,, 0) j , (24) 

where M' denotes that the site is excluded in the product. The initial condition of Eq. (122!) 
is P{0) = 1. After Laplace transformation, eq (l22|) leads to, 

M 

l-Vc{l-c)Y,Hrj)Hr„s) , (25) 



Pn{s) = —— — ^ 



s + l/ro + Ei=iCfc(^j; L j=i 
where q{fj,s) = dtexp{—st)q{rj,t). 

q {vj, t) is calculated by the perturbation expansion of exp (j^t^ ■ Hi is taken as the pertur- 
bation term. Hq conserves the number of quenchers. On the other hand. Hi defined by Eq. 
( !2T|) changes the number of quenchers by the amount of one and only the odd powers of Hi 
contribute in the perturbation expansion. The expansion parameter is proportional to c(l — c) 



instead of a/c(1 — c) given in the definition of Hi. 

A simple expression is obtained by the Fade approximation 

M 

^ G {fj, fi, s) k (fi) a/c(1 - c 



s + 1/to + ck 

Hr,,s) = , , (26) 

k{r^)G{ry,rw,s)k{ryj) 



11 s + 1/t"o + CK 

v=l w=l ' 



where k is defined by Eq. ( I23i) and G {fi, fj, s) is the Laplace transform of, 

/ M' \ / M' 



a=i / \i=i 



In the Pade approximation, higher order Green's functions such as, G (ri, r2, r3|r*4, rs, rg, s) 
defined similarly to Eq. (!27|) are ignored. As stated before, odd powers of Hi should be 
left in the perturbation expansion. Since Hi is a reaction term as shown in Eq. ( pTl) . the 
results gives minus contribution to the perturbation expansion. If we denote the complete 
solution including higher order Green's functions by qx {rj, s), the approximate solution q {fj, s) 
obeys, q {rj, s) > qr ifj, s) > 0. By combining this inequality with Eq. fl25l) . we find that the 
approximate expression is a lower bound of P/v(t). 

By making time differentiation of Eq. fl27|) . we obtain, 

d 1 

—G{fi,rj,t) = G {ri,rj,t) + CG {ri,rj,t) 

at tq 

- (1 - 2c)k (f,) G (f,, f„ t) - ckG (fi, r,-, t) , (28) 

where the initial condition is given by G {fi, fj, 0) = 6ij. C represents the operator describing 
hopping transitions, 

2d 

£G (r„ r„ t) = J2 ^7(2^) [g (fefc + r,, r,-, t) - G {f„ f„t)j , (29) 
fc=i 

— * 

where d is the dimensionality of hypercubic lattice, bk + ri denotes a nearest neighbor of the 
site "Ti and the sum is taken over all nearest neighbor sites. 

By introducing Eq. (1261) into Eq. (!25i) in the Laplace domain, the Laplace transform of the 
survival probability is expressed as, 

P^(^) = • (30) 

s + 1/ro + cft - c(l - c) ^ ^ A; (f^) G {r^, f^, s) k (r^) 

v=l w=l 

We can obtain the survival probability, P/v(t), by introducing the solution of Eq. ( l28i) into Eq. 
( !30|) and making the inverse Laplace transformation. The term with c(l — c) represents the 
effect of correlated diffusion, which vanishes in the dilute limit, c 0. In the opposite limit of 
c — 1, the factor c(l — c) again vanishes corresponding to the absence of diffusion since every 
site is occupied by a quencher. In both limits, the survival probability is given by. 



PN(t) = exp {—t/ro — cut) 
9 



(31) 



Eq. fl^ reproduces Eq. ([7]) obtained for the static quenching when all sites are occupied by 
quenchers. In the dilute limit, c — > 0, Eq. ( 1311) is also consistent with the known result of static 
quenching, Eq. (Q, when the reaction rate is small. In the intermediate concentration, the 
survival probability is influenced by G (r^, rj, t) defined by the probability of finding a quencher 
at position fi at time t when it starts from fj under the prohibition of double occupancy of a 
site. 



IV. LOCALIZED REACTIONS 



When reaction takes place only at a target site r^, k {fi) = A;o5r^,rR, the Laplace transform 
of the survival probability is expressed as, 

Pn{s) = ^ . (32) 

s + 1/ro + cko - c(l - c)koG (Fr, vr, s) k^ 

In the absence of diffusion and natural decay, Eqs. (!28l) and (l32ll reproduce Eq. (fTTl) derived 
by assuming the static quenching from the beginning. 

In the presence of diffusion, Eq. (1321) represents the approximate solution which interpolates 
between solutions in low and high limits of quencher concentrations. By substituting k (f^) = 
/co5r<>,rR5 the solution of Eq. fl28l) in the Laplace space can be expressed as 

A/^ ^ \ Go{rR,fR,z) 

G{rR, tr, s) = ^ ^ -. (33) 

1 + (1 - 2c) koGo{rR, tr, z) 



Go{fi,fj, s) is the Laplace transform of the Green's function satisfying, 

CGoifi,fj,t) + Sr-,,.^5{t), (34) 



dGo {fi,fj,t) 



dt 

where £ represents the operator describing hopping transitions given by Eq. ( l29l) . Go{fR, tr, z) 
in Eq. (1331) is given in terms of the Green's function for free random walks, GoirR, -tr, s), but the 
Laplace variable is modified as a result of the excluded volume interactions among quenchers 
and expressed in terms of the initial concentration of quenchers and the reaction rate, cko, 

z = s + (1/ro) + cko. (35) 
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When quenchers can migrate on all lattice sites including the target site, the Laplace trans- 
form of the Green's function in the absence of site blocking effects and the reaction can be 
written as 

Go (r^, rn, s) = ^i^-^U{s) (36) 

s 

in terms of the Lattice Green's function U{s) defined by 

U(s) = /"■■•/" d^k / , ^ , (37) 

where '?/'(s) = r/(s + F), and the structure factor is defined by \{k) = ^ X]jf=i cos {k ■ bj/bj . b 
denotes the lattice spacing. Eq. (132!) can be rewritten as, 

Pn{s) = J , (38) 

s + (1/ro) + c — 

1/ko + (1 - c)/ (1/Go - ck 

where we use the abbreviation. Go = Go{fji, -tr, z) and z is defined by Eq.(l35ll. Eq. (I38l) is one 
of the most important results of this paper. 

Eq. fl38|) is simplified in the Smoluchowski limit which is given by /cq — oo. In order to 
obtain the limit, we rewrite Eq. (138|) as 

pNis) = J . (39) 

' ^ ^^^"""^ ^ ""l/ko + (1 - c)U{z)/ {s + (1/ro) + r + ckoil - U{z))) 
By introducing the explicit expression of U{z), we find, 

lim U{z) = 1 and lim cfco(l - U{z)) = 0, (40) 

cfeo^oo cfco— >oo 

for any spatial dimension. In the limit of fco ^ oo (hopping-controlled limit), Eq. (l39l) is 
simplified into 

Subsequent inverse Laplace transformation yields a single exponential decay, 

PN{t) = (1 - c) exp [- (1/ro + cF) t] . (42) 
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In the limit of c = 1, the reaction site is occupied by a quencher at the initial time and the 
reaction takes place immediately in the limit of fco — > oo. The probability that the reaction 
site is not occupied by a quencher is given by 1 — c and the reaction takes place with the rate 
cT which is proportional to both the hopping rate and the quencher concentration. In the 
hopping-controlled limit, Eq. fj42l) is a lower bound of the survival probability. 

In the Smoluchowski limit of ko oo for a localized sink in 1 dimensional systems, the 
survival probability shows non-exponential decay if double occupancy of sites is allowed. The 



non-exponential 



21 



decay in the presence of site blocking effects is also predicted by some theories. 



22 



231] However, our approximate results predict the exponential decay in the 
limit of ko oo. Since our derivation involves the steepest descent approximation of Nakazato- 
Kitahara's theory and the Pade approximation of perturbation expansion, there should be a 
certain limitations on our theory. In the absence of diffusion our theory predicts the exact 
results of static quenching for localized reactions, regardless of the dimensionality of the systems. 
However, in the presence of diffusion, it gives only a lower bound of the survival probability. 
The accuracy of the approximation is worse in the limit of /cq — ^ oo in the presence of diffusion. 
The accuracy also depends on the dimensionality of the systems. We conjecture that the 
perturbation term which appeared by applying Nakazato-Kitahara's theory is large in the limit 
of /eq ^ C)0 in the presence of diffusion in one dimensional systems. 

V. SIMPLIFICATION BY ADJOINT EQUATION 



When the Green's function is not known, it is convenient to define the pair survival proba- 
bility, 

M 

f{fe,t) = J2G{ri,re,t), (43) 

1=1 

which describes the survival probability of a pair whose initial separation is given by a vector fe. 
From the equation for / (r^,t), the bulk survival probability can be obtained without knowing 
the Green's function. As shown below, the equation for / [fe, t) is simpler than that for the 
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Green's function. The initial condition is given by, 

/(f,,t = 0) = l. (44) 

/ {re,t) satisfies the time evolution equation with the diffusional operator £^ adjoint with C, 

I-/ (f,, t) = --f (f,, t) + C^f (f,, t) - (1 - 2c) k (f,) / in, t) - cKf (rl, t) . (45) 
dt To 

This is a generalization o 



the time evolution equation of the pair survival probability derived 
by Sano and Tachiya. |2J] In the absence of potential, and C are equal, £^ = C. From Eq. 
f l28|) . we obtain the following relation, 

MM MM 

dt 



1=1 j=l i=l j=l 

M M 

i=i i=i 



and after the Laplace transformation it leads to 

MM r : \ v^Af 

1 - 2c 



2^2^k{ri)G{ri,rj,s)k{rj) = ^— . (47) 

i=i j=i 

By substituting Eq. (1471) . Eq. (!30|) can be rewritten as, 

Pn{s) = ^ , (48) 

S + 1/T0+ \ -c^K + c{l-c)zJ2k{re)f{ri,s)j /(I - 2c) 

where z is given by Eq. ( l35i) . This is a generalization of the equation for the survival probability 
derived by Tachiya, by taking into account the site blocking effects. W The expression for the 

n n n 

reaction rate is known for localized reactions, which leads to, pj, 14, 15[ 

M 

7^1 (l/^o) + (1 - 2c)G'o(rj{,r^,z) 

By substituting Eq. (149!) into Eq. (148!) . we reproduce Eq. (!32|) with Eq. (!33|) . For localized 
reactions, (j(ri,r2,t) is known and the adjoint equation may not be needed. However, for 
long-range reactions, calculation of / (r*^, s) using the equilibrium initial condition can be easier 
than that of G{ffi, -tr, s) using the initial condition expressed by Kronecker's delta. 
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VI. STERN- VOLMER LAW 



In this section, we study the site blocking effects of diffusion on Stern- Volmer law. A Stern- 
Volmer plot is obtained from the fluorescence intensity at different quencher concentrations. 
The relative fluorescence intensity against tjq defined in the absence of quencher is given by, 

,11 

1 . . . 0). ,50) 

rio Jg dtPo{t) To 
By substituting Eq. fl30l) . rj/rjo is obtained as, 
T] 1 1 



(51) 



1/ro + CK- c{l- c)^^k (fv,) G (r;, r^, 0) k (f^,) 



v=l w=l 

In the Stern- Volmer plot, rjo/f] — I is plotted against the concentration, c, 

M M 

r]o/ri - 1 = CKTo - c(l - c)ro ^ (^t>) {r^-, ^w, 0) (r\„) , (52) 

11=1 UI = 1 

It increases linearly with c when quenchers are dilute. Deviation from linear concentration 
dependence of ?7o/''7 — 1 is theoretically obtained by solving the equation for G (ri, r2, s) given 
by the Laplace transform of Eq. fl2S]) . 

For the localized reactions, k (r^) = /co^r^/fl, we obtain the following equation by substituting 
s ^ hmit of Eq. ([33]) into Eq. ([52]), 



l/kQ + . 

^/Go{rR,rR,zo) - cko 

^koTo ^^^^ 



(l-c)f/(zo) 



c[l-f/(^o)] + (l/ro + r)Ao 
where s —>■ limit of z is introduced, 



zo = 1/ro + cfco- (55) 
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The Laplace variable given by Eq. flS^ includes the effect of initial concentration of quenchers 
and the reaction rate. This is a signature of the excluded volume interactions among quenchers. 
Eq. (153!) is positive since we can prove 

1 > ckoGo{rR, rn, Zq) > 0, (56) 

as shown in appendix. Eq. (1531) is an important result of this paper. 

In the static limit, Eq. fl55]) reduces to the following equation by substituting G'o(^_r, rji, zq) = 

1/^0, 

1 - c + l/(rofco) 

This is the exact result. 

In the limit of fco — ^ oo, the result in the hopping-controlled limit is obtained from Eq. (1541) 

as, 

Vo/V - 1 = (58) 
1 — c 

The above expression shows that f]o/f] — 1 increases linearly with increasing the hopping fre- 
quency, r, for any concentration. In the limit of c = 1, the target site is occupied by a quencher 
at the initial time and reaction takes place with probability 1 when ko ^ oo. In the opposite 
limit of c — ^ 0, r^o/"/] — 1 is proportional to the concentration c. By time integration, we can 
show that Eq. (!58|) is consistent with Eq. (142|) . 

For various lattices, the lattice Green's function, Eq. ( 1371) . is known. As an example, we 
consider random walks on the BCC lattice. The reaction takes place at the site rp with the 
rate ko. Without loss of generality, the target site can be taken at the origin of the lattice. 
Quenchers perform random walks on the lattice including the origin, and each site can be 
occupied at most by a single quencher. The lattice Green's function is known, 25 1 

Go {tr, tr, s) = I^Fi Q, i; 1; (r/(. + T)f^ | /(. + T).^ (59) 

When r < zo, we can approximate 2-^1 \', 1; ~ 1 + ^^/IQ, for ,^ ^ [26^] in Eq. fl59l) . and 
Eq. (j53l) is expressed as, 

vo/v - 1 = T ^T^^ • (60) 

+ 



ko 1/ro + r - (cfco/8)rv (^0 + r) 



2 
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In the case of c/cq > r/8, Eq. fl^ is further simphfied as, 

h ^ i/ro + r 

Eq. (!6T!) is the result vahd irrespective of the lattice structure since it can be derived by 
introducing the approximation, Go (tri'^^r^zq) ~ 1/(2:0 + E) which is valid when E < zq, into 
Eq. (1531) . In the reaction-controlled or static limit, we obtain Eq. (I57p . whereas, in the 
hopping-controlled limit, Eq. fl55]) is derived. Eq. flM]) interpolates between the static and 
hopping-controlled limits. 

In Fig. [H the general results of Eq. (153|) with Eq. (159!) are plotted for various values of 
and r. The simplified solutions of Eq. fl60|) overlap with those of Eq. fl53l) with Eq. fl59|) in 
Fig. 1. The further simplified solutions of Eq. (pTl) are also shown. The results of Eq. ( l6Ti) 
reproduce the general results except for the case of k^T^ = 100 and T/ko = 1 where a small 
deviation is found. The results indicate that although Eq. (ICTj) is derived under the condition 
of 8cfco > r, it is applicable in practice over a wide range. The results in the hopping-controlled 
limit of fco — > 00, Eq. fl58l) . are also shown for comparison. The results in the static limit, Eq. 
(!57l) . give the lower bound of rjo/r] for a given value of /cqTo. 

VII. DECAY KINETICS 



For localized reactions, the Laplace transform of the survival probability is obtained from 
Eq. ( l38i) . In BCC lattice the lattice Green function is given by Eq. ( l59l) . Therefore, when 
r < I/to + c/co, Eq. fl38|) is expressed as, 

Ms) = 1 , (62) 

'^~o^ 1/ko + (1 - c)/ [s + 1/ro + r - {cko/8) [z + F)'] 
where z is given by Eq. ([35D and approximation of Eq. fl5^ using 2F1 (^, |; 1; ~ 1 + .C^/16 
as ,^ — is introduced. In the case of ck^ > F/8, Eq. (1621) is simplified as, 

i^iv(5) = ^ ^ ^ . (63) 

^ + l/A;o+(l-c)/(s + l/ro + F) 
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FIG. 1: ijo/r] — 1 against concentration c. A) koTQ = 100; B) koTQ = 1; C) /cqTo = 0.01. In all panels, 
curves correspond to T/ko = 100, T/ko = 1, and T/ko = 0.01 from top to bottom. The thick solid 
lines indicate the general results in the presence of site blocking effects, Eq. ()53p with Eq. (|59p . In 
C) they overlap. A dashed line in A) indicates the approximate result of Eq. (|61|) . The other dashed 
lines are not visible, since they overlap with the solid lines. A dotted line in A) represents the results 
of hopping-controlled limit, Eq. (j58p . which is valid for kQ oo. Circles indicate the static solution, 
Eq. iH). 
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The inverse Laplace transformation of Eq. flB51) is obtained as, 

= ^M-t/ro) ^^^^ _ ^^^^ _ _ ^^^^ ^ ^g^^ 

s+ — s_ 

where 



r + fco ± v (r + - ^cVko 

s± = ^—^ • (65) 



Eq. (1641) together with Eq. (1651) is the resuh independent of the lattice structures. Eqs. ( jMl) and 
(165|) are derived under the condition cfco > F/S. Accordingly, the accuracy of the approximation 
decreases by decreasing the quencher concentration. The result in the hopping-controUed limit 
of /cq ^ oo reproduces Eq. (H2|) . 

In the reaction-controlled limit, Eq. (|62|) reduces to 



Pn{s) = J ^ , (66) 

' + ^+l/A:o + (l-c)/(s + l/ro) 
and its inverse Laplace transform is given by Eq. ffTOj) when I/tq = 0. 

For comparison, we present the conventional solution of the survival probability for target 



problem when site blocking effects among quenchers is completely ignored, 

ft 



PNit) = exp (^-c^ dtikc^(ti)^ 



(67) 



where the Laplace transform of k{t) is obtained from. 



k 

skAs) = ^ • (68) 

1 + koGo (ffi, tr, s) 

In Fig. [2], the numerically obtained inverse Laplace transform of Eq. fl38|) with Eq. fl59l) 
is compared with the conventional solution, Eqs. (|67ll and (!68l) . Excluded volume interaction 
is considered in Eq. ( l38l) with Eq. ( 159|) . whereas it is ignored and each quencher is assumed 
to migrate independently in the conventional solution. In all cases, our results indicate that 
the survival probability in the presence of site blocking effects decays faster than that of the 
conventional solution where the excluded volume interaction is absent. 

In the dilute limit, c -C 1, the difference between them is small regardless of the values 
of F/fco- When the decay is mainly controlled by reaction, namely, F/Zcq ^ 1, the result of 
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FIG. 2: The survival probability against normalized time, kot. I/tq = 0. A) T/ko = 100; B) T/ko = 1; 

C) T/ko = 0.01. In all panels, curves correspond to c = 0.1, c = 0.5 and c = 0.9 from right to left. The 

solid lines are obtained by the inverse Laplace transform of the general expression in the presence of 

site blocking effects, Eq. psp with Eq. (|59p . Dashed lines are obtained by the conventional expression, 

Eqs. (j67p and (j68p . In A) the dashed lines are invisible because they overlap with the solid lines. 

Dotted lines in C) represent the solution in the hopping-controlled limit given by Eq. (j42p . Circles 

indicate the solution in the static limit given by Eq. (jlOp . 
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the conventional solution is close to that in the presence of site blocking effects even at high 
concentration of c = 0.9. In this case, the excluded volume interaction is not important. As 
r/fco decreases, the deviation of the conventional solution from the results in the presence of 
site blocking effects increases at high concentrations. 

The initial decay of the solution of Eq. (!38|) with Eq. (159!) follows the results of static 
quenching from the uniform distribution, Eq. flTU]) . over a longer period than that of the 
conventional solution given by Eqs. (157|) and fl68p in the cases of B) and C). In the initial 
time range, the decay of the survival probability takes place from the configuration where a 
quencher is initially located in the vicinity of the target. In the presence of site blocking effects, 
the migration of quenchers is suppressed and the initial decay follows the results of static 
quenching over a longer period than that derived under the assumption of free migration in the 
absence of site blocking. It should also be noticed that the survival probability obtained by 
assuming static quenching is the upper bound of that in the presence of diffusion, whereas the 
solution of Eq. (1381) with Eq. ( l59l) is the lower bound. The exact solution should lie between 
them. 

When r/fco ^ 1, the initial time regime is approximated by the static quenching, and is 
followed by the hopping- controlled regime approximately described by Eq. ( l42l) as shown in 
Fig. El 

VIII. CONCLUSIONS 

We have investigated the target reaction problem in the presence of site blocking effects 
among quenchers. Quenchers migrate on any lattice sites until reaction takes place. Reaction 
rate depends on the distance between the quencher and the target. In the case of localized 
reactions, reaction takes place when a quencher comes to the target site. Once reaction occurs, 
the system becomes inert. 

The probability of reaction event is high if the excited target is initially surrounded by 
quenchers in close vicinity. As time proceeds, quencher configurations in which the quencher 
concentration near the excited target is low is more likely to survive than other configurations 
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in the ensemble. 

In the conventional theory, excluded volume interactions among quenchers are ignored, i.e., 
quenchers are regarded as independent of one another. However, quencher migration is hindered 
by the presence of other quenchers; the migration of a quencher is influenced by the time- 
dependent positions of other quenchers due to the site blocking effects. We take into account 
the excluded volume interactions among quenchers by applying Nakazato-Kitahara's theory 
of vacancy-assisted diffusion. Our analytical solutions of the survival probability interpolate 
between those in two limits of low and high quencher concentrations and the approximation 
is good when the intrinsic reaction rate is low. When the intrinsic reaction rate is high and 
the condition for truncation of perturbation expansion, c(l — c)kl/{ckQ + F)^ < 1, is not 
satisfied, the higher order terms in the expansion is only partly taken into account by the Fade 
approximation. In other words, the higher order correlations originating from the diffusional 
collisions are not fully accounted for in the Fade approximation. In the presence of diffusion, 
our result is a lower bound of the survival probability. In the static limit, the exact results are 
reproduced from the Fade approximation. 

The decay of the survival probability has been investigated for the target problem, where 
the target is excited by a pulse initially. The initial decay is well approximated by the static 
quenching. In particular, when the hopping frequency satisfies the relation, T/ko < 1, the 
initial decay of the survival probability at high quencher concentrations follows that of static 
quenching over a long period. The long time behavior of the general solution is approximated 
by the result in the hopping-controUed limit, Eq. fH2|) . when F/Zcq 1. 

When the intrinsic reaction rate satisfies the relation, F/feg ^ 1, the conventional results 
in the absence of the site blocking effects reproduces those in the presence of the site blocking 
effects even at high quencher concentration of c = 0.9. In the reaction-controlled limit, the 
excluded volume interaction among quenchers is not so important as that in the hopping- 
controlled limit. 

According to the Brownian dynamic simulation, the survival probability in the presence of 



excluded volume interaction decays faster than that in its absence. 
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28l | Similar enhancement 



of deactivation was also found by numerical simulation of random walk model on 1 and 2 



21 



dimensional lattices. |20|, |2l|, |22|, |29|] Our results also suggest that the decay of the survival 
probability is accelerated by site blocking effects in the case of lattice random walk. The result 
can be understood as follows. Assume that there are N quenchers on the lattice. As time 
proceeds, quenchers hop from site to site. Therefore, individual sites occupied by quenchers 
change with time. As long as the target site remains outside the sites occupied by quenchers, 
reaction does not occur. Once the sites occupied by quenchers include the target site, reaction 
occurs. In the absence of site blocking effects, different quenchers are allowed to occupy the 
same site. Therefore, in this case the number of the sites occupied by quenchers is generally less 
than A^. In the presence of site blocking effects, different quenchers are not allowed to occupy 
the same site. Therefore, in this case the number of the sites occupied by quenchers is A^. In 
other words, the number of the sites occupied by quenchers is generally larger in the presence 
of site blocking effects than in its absence at any time. Accordingly, the probability that the 
sites occupied by quenchers will include the target site is higher in the presence of site blocking 
effects than in its absence at any time. As a result, the survival probability of the target decays 
faster in the presence of site blocking effects than in its absence. 

It is interesting to note the quite opposite effect of site blocking on the survival probability 
of a geminate pair with a large initial separation. Recently, it has been shown that the pair 
survival probability decays slower in the presence of site blocking effects by inert particles, {s^ 
Here, the diffusion toward the target is just hindered by inert gases. 

Finally, we comment on the excluded volume interaction between the target and a quencher. 
It is possible to exclude the origin occupied by the target for the random walk of quenchers by 
modifying the lattice Green's function of periodic lattice. The research in this direction is now 
undertaken. 
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APPENDIX A: PROOF OF 1 > cA;oGo(rR, tr, z) > 

We first prove 1/G'o(^r, '>^r, z) — cko > which can be transformed into 1 > cA;oGo(rR, tr, z). 
By introducing 

z 

we obtain 

l/Go(ri?, tr, z) - cko = ^^^-^ . (Al) 

From the definition of U{s) given by Eq. fl37j) . we can show U{z) > 1 since the denominator 
in the integrand of U{s), 1 — tp{z)X{k), is smaller than 1. Since the denominator of Eq. ( lAip 
is positive, we need to prove the positivity of the numerator, cko (1 — U{z)) + (I/tq) + F > 0. 
By using 

1 - Uiz) = v^J--- r d'k "^^^^L.n > (A2) 



z + T 



1 - \{k) 



T + cko{l-U{z)) = ^J ... I d'k ■ \ (A3) 



(27r 

the numerator of Eq. ( lAll) can be rewritten as, 

1 f r Ji^ + V^o) + r (cfco + F) f 1 - \{k) 

s + 1/ro + cko + T(l- A(fc)) 

Since X{k) = ^ Xljli "^o^ (^^ ' ^j/^) < d^S]) is positive. Therefore, ckQ{l — U{z)) + 

(I/tq) + F > and it leads to l/GQ{fR,rR, z) — ck^ > 0, which can be rewritten as, 

1> ckoGo{fR,fR,z) > 0, (A4) 

where we have used the fact that both cko and G'o(rR, ^r, z) are positive, z in Eq. (]A4I) is given 
by2 = s + l/ro + cfco- Therefore, if we take the limit of s — in Eq. f lA4l) we have Eq. ( l56l) . 
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